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ABSTRACT: Exergy is an index to measure energy quality, 
which reveals the essence of work capacity loss in the process 
of energy transfer. In order to study the dynamic evolution 
mechanism of exergy in time-varying energy networks, the 
constitutive relations between physical quantities in energy 
networks are listed based on the energy network theory, and 
the generalized expressions of exergy are given, including the 
generalized description and loss equation of exergy based on 
the second law of thermodynamics, and the energy level 
factor is introduced to evaluate the energy quality. The 
dynamic evolution process of different forms of exergy in the 
energy transfer tube (line) is analyzed, including electric 
exergy, thermal exergy and pressure exergy, and the output 
equation and efficiency calculation method of exergy are 
given. The dynamic evolution process of exergy in energy 
conversion equipment is analyzed, and the loss, storage and 
efficiency of exergy are calculated. Finally, the dynamic 
evolution process of exergy in an integrated energy system is 
simulated by a specific example. The study of this paper can 
fully tap the energy efficiency potential of the integrated 
energy system, and lay a solid theoretical foundation for 


better realization of energy cascade utilization. 
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摘要 : 烟 是 衡 品质 的 指标 ， 深 刻 揭 万 
中 做 功能 力 损失 的 实质 。 为 研究 时 变 能 量 网 络 中 烟 的 动态 
演化 机 理 ， 本 文 基于 能 量 网 络 理论 列 写 能 量 网 络 中 物理 
间 的 本 构 关 系 ， 给 出 烟 的 普遍 化 表达 式 ， 包 括 基 于 热力 学 


:能 量 传递 过 寸 程 


ji 
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的 普遍 化 描述 与 烟 损 方程 ， 同 扩 
; 分 析 了 不 同形 式 ) 
和 态 演 化 过 程 ， 包 括 电 烤 、 热 ， 
方程 与 效率 计算 方法 ， 分 析 了 能 量 转换 
的 损耗 量 、 
人体 算 例 对 烟 在 某 综合 能 源 系 
进行 动态 仿真 。 本 文 所 研究 的 
能 够 充分 挖掘 综合 能 源 系统 的 能 效 潜 力 ， 为 更 好 地 实现 
能 源 梯级 利用 打下 了 坚实 的 理论 基础 。 
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关键 词 : ) 
系统 ， 能 效 潜力 
0 引言 

为 应 对 全 球 能 源 与 环境 危机 ， 


推进 低 碳 综合 能 源 


网 络 理论 ， 综合 能 源 


各 个 国家 积极 
源 系 统 ( integrated energy 
systems，IES ) 的 建立 ， 以 促进 能 量 梯级 利用 ， 
实现 节能 减 排 中。 因此 ， 号 需 寻 求 提 升 IES 能 效 
的 相关 技术 。 由 热力 学 第 二 定律 可 知 ， 提 升 I[ES 
能 效 不 仅 要 关注 能 量 “ 量 ”的 变化 ， 更 要 关注 能 
量 “ 质 ”的 变化 ， 即 烟 (exergy ) PE 。 研 究 籼 
的 动态 演化 机 理 能 够 揭示 能 量 传递 、 转 化 过 程 中 
有 用 能 的 存储 与 耗 散 情况 ， 发 现 IES 用 能 的 薄弱 
环节 ， 从 而 给 出 合理 的 [ES 能 效 提升 指导 意见 


可 尽 
3] 


o 


自 IES 概念 提出 ， IES 的 优化 运行 与 控制 等 
领域 已 有 诸多 研究 ,例如 文献 [9] 提出 基于 优化 
的 住宅 型 能 源 枢纽 的 能 量 管理 方案 ， 考 虑 了 热 ! 
联 产 单元 、 插 电 式 混合 动力 汽车 、 太 阳 能 电池 板 、 
鞭 热 等 部 分 以 及 电 - 气 输入 间 的 相互 影响 ， 文 献 


GD 
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[10] 对 冷 热电 联 供 的 IES 建立 经 济 性 数学 模型 ， 用 物理 量 间 本 构 关 系 及 统一 符号 ;基于 料 的 物理 
并 基于 模糊 层次 分 析 法 研究 了 系统 多 目标 优化 评 意义 ， 给 出 烟 的 普遍 化 描述 ， 结 合 拉 格 朗 日 方法 
佑 方法 与 运行 控制 策略 ; 文献 [11] 为 解决 集成 推导 出 烟 损 的 普遍 化 表达 方程 ， 引 入 能 级 因子 衡 
IES 的 经 济 性 、 排 放 性 和 多 目标 调度 问题 ， 提 出 量 能 量 品质 。 


了 一 种 结合 改进 的 交叉 粒子 群 算法 和 代替 价值 交 2 ) 分 析 其 体形 式 烟 的 动态 演化 机 理 ， 基 汪 

换 法 来 寻找 最 优 的 非 支配 解 。 然 而 ， 以 往 的 研究 时 变 能 量 传递 模型 ， 分 析 电 、 热 、 压 三 种 类 型 能 

中 大 多 只 考虑 能 量 守恒 定律 而 忽略 了 能 量 在 不 可 量 在 时 变 网 络 中 的 耗 散 与 存储 情况 以 及 能 量 中 不 

逆 传 递 过程 中 的 降 质 现象 。 可 用 能 的 大 小 ， 最 后 给 出 了 烟 的 输出 方程 和 效率 
随 着 IES 能 效 提升 技术 遇 到 瓶 诺 ， 有 部 分 学 I 

者 开始 重视 烟 在 IES 中 的 应 用 。 己 有 学 者 将 烟 指 分 析 能 量 转换 设备 中 籼 的 动态 演化 机 理 : 

标 纳入 IES 的 综合 能 效 评价 ， 文 献 [12] 引入 了 能 a 


值 理论 ， 建 立 基于 能 值 分 析 的 IES 价值 定量 评估 0 
方法 ， 有 助 于 挖掘 IES 能 源 互补 价值 ， 也 有 学 者 菜 的 损耗 、 存 储 与 做 功效 率 情况 分 析 。 


一 


将 烟 效 率 引 入 了 IES 的 多 目标 优化 调度 中 中， 建 4 ) 通过 具体 算 例 ， 对 烟 在 网 络 中 的 动态 演 
立 并 求解 了 综合 考虑 能 源 利用 效率 、 经 济 效益 以 ” ”化 过 程 进行 仿真 与 分 析 。 

及 环境 效益 的 IES 优化 调度 模型 :文献 [14] 则 给 。 1 ”能 量 网 络 中 烟 的 普遍 化 表达 

出 了 考虑 烟 的 区 域 IES 多 目标 规划 方法 ， 利 用 能 Pr ee 

源 集线器 模型 构建 并 求解 了 考虑 烟 效 率 和 经 济 性 烟 指 能 量 在 环境 条 件 下 可 转化 为 有 用 功 的 部 
的 双 层 优化 结构 。 然 而 ， 尚 未 有 文献 分 析 烟 在 能 ” ” 分， 体现 了 能 量 的 品质 。 由 热力 学 第 一 定律 知 ， 
量 网 络 中 的 动态 演化 过 程 ， 无 法 深刻 解析 IES 中 ”能量 是 守恒 的 ， 在 数值 上 大 于 或 等 于 烟 ， 热 力学 
能 量 品质 的 逐 级 递减 现象 。 第 二 定律 则 指出 ， 能 量 的 具体 转换 过 程 有 明显 的 


分 析 籼 在 IES 系统 中 的 动态 演化 过 程 需 要 精 方向 性 ， 如 电能 te ee 


确 的 IES 模型 。 目 前 ， 国 内 外 学 者 已 就 IES 建 模 。 无 法 全 部 转换 为 电能 ““。 可见， 部 分 能 量 在 传 
展开 了 诸多 而 究 。 文献 [15] 指出 键 合 图 基于 能 量 递 、 转换 过 程 中 会 不 可 避免 地 退化 为 不 可 用 能 ， 


CE 


合 已 上 


守恒 原理 能 够 对 复杂 非 线 性 的 IES 进行 统一 建 模 ， 且 过 程 不 可 逆 ， 而 这 部 分 损失 的 能 量 就 是 烂 损 。 
该 建 模 方法 虽然 结构 简明 规范 ， 但 无 法 揭示 能 量 1.1 


三 ”传递 的 动力 学 机 理 。 文献 [16] 回 顾 了 适用 于 IES ”能 量 网 络 常用 物理 量 间 本 构 关系 及 统一 符号 
〇 综合 管理 的 “能 源 枢 纽 ”概念 与 模型 ， 分 析 不 足 能 量 网 络 的 物理 量 包括 四 大 类 ，/ 义 位 移 


并 展望 未 来 发 展 需求 ， 但 该 模型 内 部 是 黑箱 模型 ， 基本 三 延 量 ) @ 、 上 广义 流 变量 《基本 广 延 量 
无 法 对 IES 进行 统一 建 模 ， 在 模型 上 仍然 相对 分 。 ”六 量 ) 9 、 广义 强度 量 和 广义 动量 X 。 它 
立 。 近 年 来 ， 借 鉴 电网 络 理论 的 IES 统一 建 模 方  ” 们 之 间 的 本 构 关系 为 : X=Rp ，@D=Cx 

法 逐渐 被 关注 加， 这 种 方法 根据 不 同 能 源 网 络 。 ”习习 ?9 。 不 同 能 量 网 络 的 物理 量 及 符号 详 见 表 
间 物 理 量 的 本 构 关系 ， 利 用 广义 阻 性 RR 、 容 性 ”【， 包 括 电力 、 热 力 、 水 力 和 机 械 能 量 网 络 中 上 
C ， 和 感性 元 件 壹 “， 建 立 通用 的 能 量 传递 模型 述 四 类 广义 变量 以 及 能 量 功 率 、 火 用 功率 以 及 能 


Pen， 模型 统一 、 简 明 ， 这 将 给 烟 的 动态 演化 分 析 。 ” 级 系数 。 


带 来 极 大 便利 。 因 此 ， 本 文 利用 基于 能 量 网 络 理 12 

论 中 建立 的 IES 统一 时 变 模型 ， 首 次 提出 一 种 时 。 ” 舌 的 普遍 化 描述 

变 能 量 网 络 中 烟 的 动态 演化 机 理 分 析 方 法 ， 精 确 gi 

度量 了 不 同形 式 能 量 在 IES 传递 过 程 中 对 应 形式 量 与 三 延 量 ， 其 中 三 延 量 又 包括 基本 三 延 量 、 能 

谢 的 存储 量 、 损 耗 量 以 及 效率 ， 具 体 研究 分 为 4 量 与 烟 。 体 系 中 任何 形式 烟 的 函数 都 可 用 强度 量 

个 步骤: 变量 与 广 延 量变 量 表示 3。 体系 的 烟 源 于 体系 与 
1 ) 给 出 烟 的 普遍 化 表达 : 给 出 能 量 网 络 常 环境 间 的 强度 量 差 ， 由 此 得 到 任意 形式 烟 的 普遍 


第 00 期 陈 思 敏 ， 陈 掉 勇 ， 周 孝 信 等 ， 时 变 能 量 网 络 中 烟 的 动态 演化 机 理 分 析 与 仿真 3 


化 微分 式 如 下 。 Xiodg@ 项 表征 能 量 中 不 能 被 利用 的 部 分 。 基 于 式 
dE, =(x; — Xio)d® ( 1) ( 1 ) ， 探 讨 能 量 微观 传递 过 程 品 质 的 宏观 变化 ， 
口 


式 中 ， 下 标 表示 与 第 i 种 形式 能 量 有 关 的 以 热力 系统 能 量 品质 分 析 为 例 ， 如 图 1 所 示 。 

量 ， 有 为 师 ，Xi 为 强度 量 的 疲 态 值 。 式 中 
表 1 不 同 能 量 网 络 的 物理 量 及 符号 

Table 1 Physical quantities and symbols of different energy networks 


变量 电力 变量 热力 变量 水 力 变量 机 械 变 量 
广义 位 移 八 电荷 四 /C 烂 急 /JK" 体积 By /ma 位 移 四 /Tad 
广义 流 变 量 9; 电流 pq /A 和 流 9 /JK's’! 体积 流量 pv / mas 速度 po / rad's: 
广义 强度 量 Xi 电势 XU /V 温度 xT /IK 压强 加 /Pa 转 矩 Xt /Nm 
广义 动量 XX; 人 磁 通 链 Xvy /Wb 温度 ”- 时 间 积 分 XT 压强 动量 Xp /Pas 动量 Xt /Nms 
/K's 
能 量 功 率 甩 电功率 及 /W 热 功率 媚 /IW 压 功 率 三 /W 机 械 功率 万 /W 
烟 功 率 Pei 电 烟 功率 fe / 热 炸 功 率 Bn /IW 压 烟 功 率 fxp /IW 机 械 烟 功率 Bj / 
W W 
能 级 系数 &i 电能 能 级 se 热能 能 级 sh 压 能 能 级 2p 机 械 能 能 级 5j 
1.3 
四 烟 损 普遍 化 方程 
让 可 结合 式 〈 1 ) 与 拉 格 朗 日 方法 P9， 得 到 烟 的 
NXT 2 Ea 动力 学 平衡 方程 如 下 。 
有 DFE =- 炎 JE + 籽 x; + (Yi Yio)o; 
ro | 
| > (2) 
3 式 中 ，Jss 是 烟 流 密度 ，J， 是 基本 广 延 量 
1 热力 系统 能 量 品质 分 析 示 意图 流 密度 ，o， 是 基本 广 延 量 源 强度 。 等 号 左边 为 
Fig.1 Schematic diagram of energy quality analysis of 单位 体积 内 烟 随 时 间 的 变化 量 ， 等 号 右边 第 一 项 
thermal system 为 净 流 入 体积 元 边界 的 烟 ， 第 二 项 为 在 强度 量 差 


没有 定数 量 的 热能 从 系统 1 的 恒定 热源 。 推动 下 不 同形 式 烟 间 相互 转换 后 产生 的 烟 ， 第 三 


省 x。 保持 不 变 。 由 热力 学 第 一 定律 可 知 由 于 综合 能 源 系 统 中 能 量 通 常 在 圆柱 形 
TO BE o 了 ? 


~ 


二 a | 加 sem 页 大 ) 正 量 源 产 避 ， 0 
加 传递 到 系统 II “的 恒定 热源 zm。， 温 度 寂 态 ”项 和 随 基本 三 延 量 源 产 生 的 类 本 


XriD1 = XpD。; 由 热力 学 第 二 定律 可 知 ， 随 着 eA i ee 
过 程 的 进行 精 不 断 增 大 ， 不 可 用 能 也 随 之 由 您 积 分， 并 移入 整理 后 得 到 下 3 


(XA Xo) Pa (xs Xo ) Pe= Ce pdx 
系统 I 的 烟 Bu 。 分 析 可 知 ， 传 热 过 程 发 生 烟 损 ， dt 

热能 做 功能 力 损失 了 ， 热 能 降 质 。 可 见 ， 要 保证 (x -Xo) dx 
能 量 数量 不 变 就 要 牺牲 能 量 的 品质 ， 反 之 ， 要 维 

持 能 量 品质 不 变 就 要 损失 更 多 数量 的 能 量 。 9 


式 中 ， 省 略 了 下 标 i ， 角 标 A、B 分 别 代表 
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恒定 律 知 ， 电 荷 源 强度 04=0， 由 式 (4) 得 到 


能 量 传 入 截面 A 与 能 量 传 出 截面 B，x。 是 能 量 


电 烂 功率 损失 方程 如 下 。 
传递 过 程 与 长 度 相关 的 基本 广 延 量 的 函数 ， AP。=AR + xuo (Ps -pos)=APR (6) 
方程 (3 ) 等 号 左边 两 项 即 单 位 时 间 能 量 在 

截面 A、 0 如 下 。 式 中 ，AP. 是 电能 从 截面 A 传递 到 截面 B 产 
生 的 电 燃 损 ，AP 是 传递 过 程 的 电能 损失 ， Xu 
AR =AP+ tA- (XY 一 加 ) i。 为 电势 寂 态 值 ， 由 于 通常 取 大 地 为 零 电势 点 ， 故 
(4) Xuo=0 ， 电 能 能 级 s。 =1 。 因 此 ， 电 烂 损失 与 电 

式 中 ，AP. 是 能 量 从 截面 A 传递 到 截面 B 产 。 ”能 损失 变化 相同 。 
生 的 烟 损 ，AP 是 传递 过 程 消 耗 的 可 用 能 ， 以 电能 在 平行 双 传 输 线 传递 过 程 为 例 分 析 电 


NAP 代表 了 不 可 逆 性 产生 的 可 用 能 损失 。 可 见 ， ”和 烟 的 演化 ， 时 变 等 值 模型 如 图 2 所 示 。 图 2 中 ， 
能 量 传递 过 程 产生 的 烟 损 主要 包含 了 传递 过 程 消 及 是 电阻 ， 由 于 电阻 发 热 使 得 部 分 电能 转化 成 
耗 的 可 用 能 以 及 不 可 逆 性 造成 的 可 用 能 损失 ， 又 。 ” 低 品 位 热能 ， 产 生 烟 损 ;，C。 是 电容 ， 大 小 与 电 
三 等 于 体积 内 烟 随 时 间 的 变化 量 减 去 强度 量 差 推动  ” 势 X0 变化 快 侯 有 关 ， 体 现 烟 的 存储 效应 ， 工 . 
站 下 转化 为 该 形式 的 烟 ， 再 减 去 随 基本 广 延 量 源 而 


是 电感 ， 反 应 电流 9。 变化 的 快慢 ， 体 现 烟 的 惯 
产生 的 籼 。 
1.4 性 效应 。 电 烟 功 率 损失 计算 方法 与 电能 功率 相同 ， 
能 级 根据 电路 理论 外 可 以 求 得 ， 此 处 不 更 述 。 
引入 能 级 因子 评价 能 量 的 品质 。 体 系 中 第 
i 种 形式 能 量 对 应 的 能 级 因子 s。 的 大 小 等 于 对 A 0 
应 烟 微 分 与 能 量 微分 ( dE, =xid@ ) 的 比值 ， es c, 2 
si 大 小 如 下 。 
2i = 二 ($5) 
t i 2 电能 传递 过 程 的 时 变 等 值 模型 
分 析 可 知 ， 当 ， 表 明 此 时 体系 强度 量 Fig.2 Time-varying equivalent model of power transfer 
Xi ee 
三 。 ”表明 强度 量 寂 态 值 Z。 为 0， 能 量 在 数值 上 等 于 。 ，。 
© 烂 ， 此 时 体系 烟 达 到 最 大 值 。 可 见 ， 体 系 能 级 越 


ee 、 热 烟 的 动态 演化 
高 表明 能 量 中 烂 的 比例 越 大 ， 能 量 作 用 于 相关 外 热力 网 络 中 管道 流体 在 压力 差 的 推动 下 发 后 


界 的 功 就 越 多 ， 则 能 量 的 利用 价值 越 显著 。 对 流传 热 。 热 烟 表 征 由 于 体系 温度 与 环境 温度 的 

2 ”具体 形式 烟 的 动态 演化 机 理 不 同 而 能 向 相关 外 界 作 最 大 有 用 功 的 能 力 。 由 式 
按照 能 量 网 络 中 物质 的 运动 形式 分 类 ， 得 到 “4) 得 到 热 话 功率 损失 方程 如 下 。 

不 同形 式 的 烟 ， 例 如 ， 电 力 网 络 中 电荷 定向 移动 

产生 电 烂 ， 热 力 网 络 中 业 的 运动 产生 热 烟 ， 水 力 

网 络 中 体积 流量 的 运动 产生 压 烟 。 以 下 基于 能 量 


网 络 理论 所 建 时 变 模型 ， 分 析 具 体形 式 烟 在 管 dB, 8 
APr =AR + XroAP.,= — 9.dxr 一 
( 线 ) 控制 体 截面 A、 截 面 B 间 的 动态 演化 机 理 。 dt A ( 7) 
2.1 (Kr — Xro ) ds 
电 烟 的 动态 演化 l 
电力 网 络 中 ， 电 燃 表 征 电场 作用 下 带电 体系 式 中 ，APs 是 热 蔽 损失 ，AR 是 热能 损失 ， 


向 相关 外 界 做 出 的 最 大 有 用 功 的 能 力 。 由 电荷 守 9, 是 人 流 ，Xs 是 炉 源 随 管道 长 度 变化 的 函数 ， 


2 A 
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时 变 能 


三 


量 网 络 


日 ， 


的 动态 演化 机 理 分 析 与 仿真 5 


Mro 为 温度 彼 态 值 ， 通 常 取 为 环境 温度 ， 
(5) 可 知 热能 能 级 sr 小 于 1， 热 能 属于 低 品 


由 式 


立 能 


子 


由 式 〈7 ) 分 析 得 ， 热 烟 损 失主 要 包括 热能 
损失 以 及 传 热 和 扩散 不 可 逆 性 引起 的 能 量 损失 ， 
也 等 于 单位 时 间 内 总 热 灶 的 变化 量 〈 非 稳 态 项 ) 
减 去 热 烤 的 产生 量 。 以 一 维 单 层 圆柱 形 管道 中 热 
能 的 传递 过 程 为 例 具 体 分 析 。 由 能 量 网 络 理论 建 
模 方法 得 到 传 热 过 程 等 值 回 转 器 GY 模型 如 图 3 
所 示 。 


图 3 对 流传 热 过 程 的 时 变 等 值 模型 
Fig. 3 Time-varying equivalent model of convective heat 
transfer process 
中 ，Zxra 是 输入 温度 ， xirs 是 微 元 控制 体 
内 的 温度 ，@ 是 热 容 ，R 是 管 壁 径 向 总 热 阻 。 由 
图 3 所 示 模 型 ， 写 出 传 热 微 元 控制 体内 的 能 量 损 


耗 如 下 。 
AR =C Mrp (a) 
AP ~ (Xm Xm) (b) 


AR,=(1 — Er ) Ps -(1-— éna ) Ba (c) 


= XTo (pa 一 DA) 


(8) 

式 中 ，AR. 是 储 热效应 下 热 容 元 件 中 存储 的 
热 功 率 ， 即 流体 温度 随时 间 变 化 而 在 管道 累积 的 
能 量 ，AR, 是 热 阻 元 件 上 损失 的 热 功 率 ， 即 导热 
与 对 流 过 程 沿 管 侧 泄漏 的 热能 ，AR。 是 热能 
日 于 能 级 下 降 所 产生 的 不 可 逆 耗 散 对 应 的 热 功 率 
损失 ，shA=1-Mro/NrA 是 控制 体能 量 输入 处 热 
能 能 级 ，éhs=1- Xro/Xrs ”是 控制 体内 热能 能 级 。 

结合 公式 (7 ) 、〈 8 ) ， 得 到 控制 体内 热 
烟 的 动态 演化 方程 如 式 ( 9 ) 所 示 。 


Pus—Paa 


全 


EN Wi Xro) 
h 


Ch XrB XTo (pa 一 DA ) 


(9) 
式 中 ， Prs 为 控制 体 输出 的 热 烤 功 率 ， Rs 
为 控制 体 输入 的 热 烟 功率 。 由 式 〈9 ) 画 出 热 烟 
在 单 层 圆柱 形 管道 传递 过 程 的 动态 过 程 如 图 4 所 
示 。 


XTO (9sB = OAA) 


号 xhA FpB 


Ch 多 TB 


(xre — Xro0)/Rr 


图 4 单 层 圆 柱 形 管 道 热 烟 动 态 演化 示意 图 
Fig. 4 Diagram of thermal dynamic evolution of single- 


layer cylindrical pipe 


4 易 得 ， 热 烟 在 管道 中 传递 的 烟 效 率 如 


由 图 


(Ma 一 Xro) 
A 
Fra 


Ch XrB KTo (pm 一 pwA) 


(10) 
2.3 
压 烟 的 动态 演化 
在 流体 网 络 中 ， 压 烂 表征 体系 在 不 同 压 力 环 
境 下 对 相关 外 界 作 最 大 有 用 功 的 能 力 。 流 体 网 络 
中 体积 源 强 度 cv =0， 由 式 〈4) 得 到 压 烟 功率 
损失 方程 如 下 。 


AP,=AR + XpoAP, = 


B 
— Pvdxp 


dt 

(11) 

式 中 ，AAPu 是 压 能 传递 过 程 压 烟 功 率 损 失 ， 

A 下 是 压 能 传递 过 程 的 压 能 功率 损失 ，9, 是 体积 

流量 ， Xe 为 压强 披 态 值 。 

以 圆 形 管 路 的 流体 压 能 传递 过 程 为 例 具 

析 压 籼 的 演化 ， 由 能 量 网 络 理论 建 模 方法 得 到 压 

能 传递 过 程 通用 等 值 模型 如 图 5 所 示 。 考 虑 在 一 

维 管 流 中 ， 平 铺 的 流体 管线 ， 法 向 方向 上 流体 受 

到 的 压力 和 重力 相 平衡 ， 管 路 和 流体 不 发 生 弹 性 
形变 与 线性 形变 ， 且 一 般 忽 略 对 流 项 。 
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5 压 能 传递 过 程 的 时 变 等 值 模型 
Fig. $ Time-varying equivalent model of pressure energy 


transfer process 


Na 


为 控制 体 输入 的 热 烂 功 率 。 从 能 量 的 角度 进 
分 析 容 性 元 件 与 感性 元 件 上 的 : 流 容 表征 流体 压 
缩 性 ， 与 压力 变化 快慢 有 关 ， 流 容 如 同 电容 将 外 
界 能 量 转变 为 电场 能 储存 起 来 一 样 ， 把 流体 的 动 
能 转变 为 压力 势能 存储 起 来 ， 当 流体 处 于 稳 态 和 
低频 时 可 不 予 考虑 ， 流感 表征 流体 惯性 ， 当 流体 
商 频 志 动 时 在 惯性 作用 下 使 流体 流量 发 生变 
速 ， 并 进一步 引起 压力 的 变化 ， 此 时 ， 流 感 如 同 
电感 将 电场 能 转化 为 磁场 能 储存 起 来 一 样 ， 把 压 


图 中 ， xps 为 截面 A 压强 ， Xe 为 截面 B 压 
强 ，wp， 为 体积 流量 ，C， 是 流 容 ，R， 是 流 
阻 ，L， 是 流感 。 由 图 5 所 示 模 型 ， 写 出 压 能 传 
递 过 程 微 元 控制 体内 的 能 量 损耗 如 下 。 


Am = XpaC, MA (a) 
AP. -9 BR, (b) 
A = 9,pL,Pv (c) ( 12 ) 


AP, =(1-sejBs-(L-sAjR (d) 


= Xpo (Pa OA) 


式 中 ，Am. 是 流窜 上 存储 的 压 功率 大 小 ， 


Ai 是 流 阻 上 产生 的 沿 程 压 功率 损失 ， 即 由 于 流 


体 的 秋 性 或 管 壁 的 粗糙 引起 的 ，Am 是 流感 存储 


的 压 功 率 大 小 ，Amm 是 能 级 下 降 的 不 可 逆 过 程 


所 产生 的 压 功 率 损 失 ，shA=1- Xeo/XNrA 是 压 能 
输入 处 能 级 ， EvB™—1— Xpo/ Xpe 是 压 能 输出 处 能 
级 。 

结合 式 (11) 、(12) ， 得 到 控制 体内 压 烟 
的 动态 演化 方程 如 下 。 
Pys = Popa AR. AP. A AP, 


= = XpACv XpA — Oo 及 -alps — (13 
Npo (pa 一 DA ) 
) 


式 中 ， hps 为 控制 体 输出 的 热 烟 功 率 ， pa 


感 
力 能 转换 为 动能 储存 起 来 。 可 见 ， 流 体 传递 过 程 
:能 与 动能 相互 变换 ， 过 程 无 损 可 逆 。 

由 上 式 得 压 烂 在 管道 中 传递 的 烂 效 率 如 下 : 


A PACyXpA 一 vB 和 
几 -i x DA ea 
小 人 — Xpo (pa x Oa ) 
3 能 量 转 换 设 备 中 类 的 动态 演化 


能 量 转换 设备 是 耦合 不 同 能 源 子 网 络 的 重要 
元 件 ， 最 为 常见 的 类 型 有 感应 电动 机 和 离心 泵 ， 
以 下 基于 能 量 网 络 理论 所 建 时 变 模型 ， 分 析 烟 在 
其 中 的 演化 机 理 。 

3.1 
感应 电动 机 

通过 电磁 感应 ， 感 应 电动 机 将 输入 电压 x 
与 感应 电流 m 共同 作用 所 产生 的 电功率 号 转换 成 
感应 转 矩 x 与 线圈 转动 速度 流 9。 共同 作用 所 产 
生 的 机 械 功 率 。 机 械 能 是 高 品质 能 量 ， 能 级 系数 
为 1， 因 此 机 械 类 功率 记 等 于 机 械 能 量 功 率 。 

根据 能 量 网 络 理论 ， 将 感应 电动 机 等 效 为 回 
旋 器 GY(k,) ， 得 到 的 时 变 等 值 模 型 如 图 6 所 示 。 


旺 


Pa 


6 感应 电动 机 时 变 等 值 模型 
Fig. 6 Time varying equivalent model of induction motor 
图 中 ， Rs 表征 电动 机 转速 差 ， Re 表征 电 
动机 轴承 与 转子 摩擦 ;Zu 表征 电动 机 转子 惯性 ; 
Cx 表征 机 轴 的 柔 度 等 特性 ， ,为 电动 机 转 矩 转 
换 系 数 ( Xi /k=%9 ，9, 使 ， Xt:， 用 于 能 量 方 
程 列 写 ) 。 由 图 6 分 析 得 到 能 量 转换 过 程 产 生 的 烟 


本 
ml 
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损失 如 下 。 | APri 0°R (a) 
Co (a) APAP Lo  (b) 
No As OO) 
APiT xo 9 (c) AR = 00R， (dj 
AP.=xCaxX: (dg 


式 中 ，Afi。 是 转动 机 械 功率 损失 ，A:i 是 摩 
擦 机械 功 率 损 失 ，AA 是 感性 存储 机 械 功 率 ， 
AFi. 是 容 性 存储 机 械 功 率 ,根据 图 6 分 析 可 得 ， 感 


应 电动 机 输出 的 机 械 籼 功率 有 方程 以 及 感应 电动 
机 的 烟 效 率 me 计算 公式 如 下 。 


Fi 一 到 入 /Re 9 多 mx@10 到 [a 让 XC 万 
(16) 


3 玉 和 /Ra 0 名 mx LA 内 ’ XC 万 
G 
P 


[3 


(17) 
3.2 
离心 泵 
离心 泵 将 转 矩 办 和 角速度 流 8。 产生 的 机 械 
功率 如 转换 成 流体 压强 加 和 流量 wm, 产生 的 流体 
压 功率 中， 使 流体 升 压 。 离 心 泵 转动 部 件 的 机 械 
动机 提供 ， 使 流体 压 能 增加 的 烟 称 


H 


| 


时 网 络 理论 ， 将 离心 泵 等 效 为 
GY (1, ) ， 得 到 离心 泵 的 等 值 


旋 器 
模型 如 图 7 所 示 。 


L: Re 9, 


上 


Nt GY Kp 


7 离心 泵 的 时 变 等 值 模型 
Fig. 7 Time varying equivalent model of centrifugal pump 
图 中 ， Ar 表征 转子 摩擦 特性 ， Lox 表征 转子 
转动 惯量 特性 ， 如 :表征 泄漏 特性 ， 心 :表征 阀门 
出 口 压力 损失 特性 ， 丈 为 叶轮 回旋 器 模 数 
/=9,，9om 会 px ， 用 于 能 量 方程 列 写 ) 。 
由 图 7 分 析 得 到 能 量 转 换 过 程 产生 的 烟 损 失 为 : 


LH 


| 


式 中 ，AP' 是 叶轮 机 械 功 率 损失 ， 主 要 由 离 
心 泵 转子 摩擦 造成 ，API 是 存储 的 机 械 功 率 大 
小 ，Ah 是 容积 功率 损失 ， 主 要 由 离心 泵 的 泄漏 
特性 造成 ，AB 是 水 力 功率 损失 ， 主 要 由 阀门 出 
口 压力 损失 造成 。 

考虑 到 液体 压强 寂 态 值 不 为 零 ， 部 分 压 能 
法 转化 为 有 用 功 ， 并 产生 不 可 着 能 量 损 失 Neopw 。 
故 机 械 能 在 离心 泵 中 转换 为 液体 压 能 时 ， 离 心 泵 


输出 的 压 烟 功率 籼 A， 方程 与 离心 泵 的 烟 效 率 
公式 如 下 。 


2 
PL 一 Ps? RA, 


xp 


可 


mB 计生 
— 1/ Bovey RR -pL 9 -X09 
(19) 


_ Ps? R, — Xl/ Rovete AR 一 人 L 9 = 9 
Pp 


] 


(20 ) 

由 式 (16) 与 (19) 得 到 IES 能 量 转换 元 件 

将 电能 转化 为 压 能 的 过 程 中 燃 的 动态 演化 如 图 8 
所 示 。 


图 8 感应 电动 机 与 离心 泵 中 烟 的 动态 演化 示意 图 


Fig. 8 Dynamic evolution diagram of exergy in induction 


motor and centrifugal pump 


4 ” 算 例 分 析 络 。 
4.1 e; eI ©, e4 
算 例 描述 1; 
本 节 以 某 综合 能 源 系统 为 例外 ， 拓 扑 图 如 图 hp3 hl hp hl np 
9 所 示 ， 该 系统 含 4 节点 电力 网 络 以 及 6 节点 热 Re 
力 网 络 ， 由 热电 联 产 机 组 (CHP) 分 别提 供电 源 
与 热源 ， 热 电 比 为 1.3 。 电 力 网 络 节 点 6 连接 (a) ”电力 网 络 和 热力 / 水 力 网 络 
电源 ， 节 点 e 、e) 负荷 大 小 都 为 0.15 MW， (a) Power network and thermal /hydraulic network 


节点 6 为 连接 大 电网 的 松弛 节点 。 热 力 供水 网 


63 ni n, ns; ns n 
络 节点 h， 、h， 负荷 大 小 为 0.3 MW ， 节 点 h， Et 
连接 热源 ， 回 水 网 络 对 称 。 供 、 回 水 网 络 间 连 接 . 
热源 与 循环 水 泵 ， 循 环 水 泵 由 感应 电动 机 驱动 。 ” " 

©o no ne 


热电 联 产 机 组 63 el e2 e4 


(b) 能 量 耦 合 元 件 


电网 (b) Energy coupling element 
= 10 某 区 域 综合 能 源 系统 能 量 网 络 的 图 
S mi! ph |/ dh2 Fig. 10 Diagram of energy network of a regional 
es 。 a | 风 comprehensive energy System 
Cy 热电 联 产 机 组 Xs rol To 加 et 
WO ,uf . 在 流体 网 络 中 ， 根 据 储 能 元 件 个 数 选取 状态 
CN 变量 : 选取 热力 网 络 管道 末端 温度 ， 水 力 网 络 管 
QO ss ~ 道 流量 与 管道 末端 压强 。 共 16 个 状态 变量 (包含 
本 和 5 个 温度 变量 ， 4 个 压强 变量 ， 6 个 体积 流量 变 
人 量 ) ， 在 网 络 中 可 以 列 出 16 个 方程 。 结 合 能 量 耦 
< Fig. 9 Topological diagram of a regional comprehensive 合 元 件 ， 可 以 再 列 出 3 个 方程 ， es 
energy system 量 ， 19 个 方程 。 在 电力 网 络 中 设 节点 、e) 
算 例 中 流体 介质 为 工 质 水 ， 电 势 秽 态 值 yo 2 个 电压 幅 值 变量 ， 除 节点 ee 
为 0Y， 压 强 彼 态 值 zz。 取 0.1 MPa ， 温 度 彼 态 。 ”变量 。 结 合 广义 基 尔 霍 夫 定律 与 烟 功 率 方程 得 到 


看 


设 为 环境 温度 xio 为 10"C 热 负 荷 出 口 处 温度 人 Re 
io 为 50 ‘C， 热 源 供应 温度 如。 为 99.9 'C。 采 用 状态 方程 () 万 () 与 输出 方程 8( ) 如 下 。 


分 段 集中 参数 模型 。 控 制 CHP 热 源 出 口 温 度 不 i 

变 ， 步 长 0.01s ， 仿 真 时 长 200 s 。 设 热 负荷 P= /pp) 

由， 由 0.3 MW 突变 到 0.32 MW ， 基 于 Matlab 软 X= 户 (PN FE 

件 ， 首 先 构建 时 变 能 量 网 络 方程 ， 利 用 四 阶 龙 格 R= 8(0,0,X, Yt) 

库 塔 法 〈 Runge-Kutta ) 求解 状态 变量 ， 结 合 4.3 

文 所 提 烟 分 析 方 法 进行 仿真 与 讨论 。 仿真 与 分 析 

4.2 1 ) 首先 给 出 流体 网 络 状 态 变 量 动态 响应 

能 量 网 络 方程 建立 te aren a 
借鉴 电路 图 论 知 识 ， 由 图 9 画 出 对 应 的 能 量 M_out 、 流 入 热 荷 2 质量 流量 M in 与 热源 返 

网 络 的 图 如 图 10 所 示 。 热 力 网 络 与 水 力 网 络 的 图 温度 Tr、 热 负 蓓 2 供应 温度 Ts load 。 


是 一 样 的 (共用 节点 与 支 路 ) ， 但 所 包含 的 内 容 
(元 件 ) 是 不 同 的 ， 图 10(a ) 中 虚线 表示 回 水 网 
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Sl | ER 一 未 源 供应 孜 疯 PXI 
班 2.9f 三 一 一 热 负荷 2? 消耗 热 烟 Pxh_load 
ee z0 
坚 1.5| 吊 
VIP 
140 Ey 而 页 120 
时 间 (s) 1 
97.3 0 30 60 90 120 
49.15 时 间 (s) 
A 49.14 
只 971 负 答 2 供应 汤 度 Ts load| | 49.13 纪 0.065 二 一 一 
二 DA 妆 一 一 热源 处 不 可 用 能 Ph0_s 
编 970 一 | 三 0.060 热 负 荷 2 处 不 可 用 能 Ph0_load| 于 
969 Sn 热 网 管道 能 量 损耗 Ph pipe 
] 49.10 呆 0. 
1 I 
2680 30 60 90 T2009 0035 
时 
对 间 (S) 
图 11 流体 网 络 状态 变量 响应 0 30 60 90 120 
时 间 (s) 
Fig. 11 Fluid network state variable response 
图 11 得 ，M_out 从 2.99 kg/s 增 加 了 3.3 % ， 图 12 热力 网 络 中 烟 的 动态 演化 图 
稳定 在 3.09 kg/s; M in 则 从 初 值 1.52 kg/s 动态 Fig. 12 Dynamic evolution diagram of exergy in thermal 


上 升 6.4 % 并 稳定 在 1.62 kg/s 。 而 Tr 由 49.125 °C 
下 降 到 49.101 "C， 再 缓慢 上 升 ， 经 过 100s 最 后 
稳定 于 49.151 "C; ”Ts load 由 97.14 'C 下降 到 
96.82 "C ， 最 后 经 过 90 s 稳定 在 了 97.18 "C ， 两 


处 温度 变化 非常 小 ， 增 量 都 不 到 0.1 % 。 热 力 网 
络 惯性 大 ， 动 态 响应 时 间 较 长 。 


分 析 得 ， 热 力 网 络 中 ， 热 源 、 热 负荷 处 温度 
变化 极 微 导致 对 应 节点 能 级 大 小 基本 保持 不 变 ， 
而 工 质 水 质量 流量 的 变化 幅度 则 是 影响 负荷 消耗 
热能 数量 与 热源 供应 热能 数量 的 主要 因素 。 

2 ) 在 热力 网 络 中 ， 关 注 负荷 突变 量 办 
处 、 热 源 处 以 及 热 网 管道 中 谓 的 演化 情况 。 热 力 
网 络 中 烂 的 动态 变化 仿真 结果 如 图 12 所 示 ， 包 括 
热源 处 供应 的 热 烤 Pxh_s 、 热 负荷 加， 消耗 的 热 
烟 Pxh_load 以 及 热源 处 的 不 可 用 能 Ph0_s 、 热 负 
荷 办 ， 处 的 不 可 用 能 Ph0 load 和 热力 网 络 管道 中 
的 能 量 损耗 Ph_pipe 情况 。 
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network 
由 图 12 可 见 ， 热 源 供应 热 烟 Pxh s 从 


0.57MW 稳定 到 0.59 MW ， 增 加 了 3.17 %; 热 负 
荷 如 ， 消 耗 热 烟 Pxh load 从 0.269 MW 稳定 到 
0.287 MW ， 增 加 了 约 6.68 %; 热源 不 可 用 能 
Ph0_s 由 0.0635 MW 稳定 在 0.0657 MW ， 增 加 了 
3.46 %; 热 负 荷 不 可 用 能 Ph0 load 由 0.030 9 
MW 稳定 于 0.032 9MW ， 动 态 上 升 6.52 % 。 管 
道 能 量 损耗 Ph pipe 由 0.034 3 MW 稳定 在 
0.0369 MW ， 增 大 了 7.58 %。 

分 析 上 述 仿真 结果 可 知 ， 热 源 处 质量 流量 
加 的 百分比 大 于 返回 温度 增加 的 百分比 ， 因 此 其 
供应 热能 增加 ， 由 于 热源 供应 温度 不 变 ， 热 能 能 
级 不 变 ， 故 热源 供应 热 烟 Pxh_s 以 及 热源 不 可 用 
能 Ph0_s 都 随 热 能 的 增加 而 增加 ;而 热 负荷 处 供 
应 温度 的 升 高 使 热能 能 级 变 大 ， 热 烟 损 失 占 比 减 
小 ， 但 其 减 小 幅度 小 于 热 负 蓓 升 高 幅度 ， 因 此 ， 
热 负荷 向， 热 籼 消 耗 量 Pxh_load 增 加 的 同时 不 可 
用 能 Ph0_load 也 增加 。 此 外 ， 源 处 的 不 可 用 能 高 
于 管道 能 量 损 耗 Ph_pipe ， 热 荷 办 ， 处 的 不 可 用 
能 接近 管道 能 量 损耗 Ph_ pipe ， 可 见 ， 能 量 网 络 
中 仅 关 注 网 络 能 量 损耗 是 不 够 的 ， 还 需要 关注 不 
可 用 能 的 大 小 。 从 本 例 仿 真 可 以 看 出 该 部 分 能 量 
占 比 高 于 管道 能 量 损耗 ， 不 应 忽视 。 

3 ) 为 使 更 多 工 质 水 产生 流动 以 满足 突变 的 
负 蓓 需求 ， 离 心 泵 出 口 处 需要 提供 更 多 的 压 籼 。 
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图 13 为 水 力 / 电力 网 络 以 及 能 量 耦 合 元 件 中 籼 的 
动态 演化 图 ， 包 括 感应 电动 机 供应 电 烟 Pxe_s 、 
离心 泵 出 口 压 烟 Pxps 、 管 道 压 能 损耗 
Pp_pipe 、 能 量 耦 合 元 件 能 量 损耗 P_coupling 以 
及 电力 网 络 能 量 损耗 Pe_pipe 。 
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于 0.00448 上 天 
双 000272 [7 二 60xl0 疆 
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13 水 力 / 电力 网 络 以 及 能 量 厢 合 元 件 中 烟 的 动态 演 
化 图 
Fig. 13 Dynamic evolution diagram of exergy in 
hydraulic / electric power network and energy coupling 


elements 


图 13 可 见 ， 离 心 泵 出 口 压 烟 Pxp_s 从 2 618W 


流 经 电阻 转化 为 低 品 位 的 热能 ， 产 生 烟 损耗 ; 
( 2 ) 热力 网 络 中 热能 含有 的 不 可 用 能 损失 到 环 
境 中 后 ， 剩 下 的 即 热 糯 ， 由 于 热 阻 与 热 容 的 存 
在 ， 在 热 网 中 随 着 流体 的 流动 发 生 损 耗 ， 到 达 负 
荷 处 又 有 一 部 分 新 的 不 可 用 能 损失 到 环境 中 ; 
( 3 ) 在 能 量 耦 合 元 件 中 ， 感 应 电动 机 消耗 电 
能 ， 在 电能 转换 为 机 械 能 的 过 程 中 伴随 着 机 械 能 
损耗 ， 感 应 电动 机 驱动 离心 汞 后 ， 机 械 能 转化 为 
压 能 的 同时 设备 产生 一 定 压 能 损耗 并 在 离心 泵 出 
处 伴随 着 不 可 用 能 的 损耗 ， 最 终 离心 倘 出 口 压 
菜 带 动 流 体 的 流动 满足 负 蓓 需求 ，〈 4 ) 在 流体 
网 络 中 由 于 流 阻 、 流 感 以 及 流 容 的 存在 产生 管道 
上 的 压 能 损耗 ， 相 比 于 热 烟 损 ， 管 道 流动 产生 的 
压 能 损 几 乎 可 以 忽略 不 计 ， 但 驱动 流体 流动 的 机 
泵 是 不 能 忽略 的 。 该 综合 能 源 系 统 用 能 过 程 中 ， 
在 能 级 基本 不 变 的 情况 下 保证 了 能 量 的 品质 ， 但 
损失 了 更 多 数量 的 能 量 ， 对 应 的 也 产生 了 更 多 的 
不 可 用 能 。 
4) 图 14 为 热力 / 电力 网 络 以 及 能 量 耦 合 元 件 
中 烟 的 总 传递 效率 。 热 力 网 络 中 热 烤 总 交 
值 为 94.4 %， 由 于 系统 储 能 元 件 的 综合 作用 ， 效 


上 升 到 2 709 WwW ， 增 加 了 3.5 % 。 离 心 泵 的 压 籼 


率 动态 瑚 荡 较 大 ， 最 后 稳定 在 94.2 %， 减 少 了 


由 感应 电动 机 电 烟 转换 而 来 ， 从 图 中 看 到 电 烦 
Pxe s 从 4527W 降 到 4478W ， 降 了 1.08 %。 
同时 ， 管 道 流 量 增加 带 来 更 多 的 管道 压 能 损耗 
Pp_pipe ， 在 动态 过 程 初始 时 刻 ， 由 于 管道 末端 
的 压强 随时 间 不 断 减 小 ， 导 臻 容 性 元 件 吸收 负 功 
率 ， 且 其 绝对 值 大 于 流 阻 与 流感 所 吸收 的 正 功 
率 ， 管 道 压 能 损耗 Pp_pipe 曲线 在 一 开始 是 小 于 
零 的 。 而 管道 压 能 损耗 从 初始 值 2.8 W 最 终 稳定 
在 3.1 W ， 基 本 可 忽略 不 计 。 此 外 ， 碍 合 元 件 的 
总 能 量 损耗 P_coupling 的 初 值 为 1586.7 W ， 并 
最 终 稳 定 在 了 1 587.2 W ， 前 后 相差 不 大 ， 但 由 
于 储 能 元 件 作 用 使 过 程 震荡 较 大 。 随 着 热源 供应 
热能 的 上 升 ， 电 源 供应 电能 也 随 之 增加 ， 电 力 网 
络 能 量 损 耗 Pe_pipe 增加 ， 从 0.033 7 MW 增 大 
到 0.036 9MW， 增加 了 9.5 % 。 

分 析 可 知 ， 热 负荷 的 增加 使 得 热力 、 水 力 、 
电力 网 络 以 及 能 量 耦 合 元 件 中 的 能 量 损耗 都 相应 
的 增加 ， 且 不 同 网 络 、 元 件 的 能 源 输入 与 输出 也 
产生 相应 变化 。 系 统 提 供 的 总 能 量 包括 CHP 机 组 
产生 的 总 电能 与 总 热能 。 (1 ) 电力 网 络 中 电能 


Tt 


0.2 % 。 热 荷 和 ， 输 出 的 热 烟 增 量 较 输入 的 热 娴 
增 量 大 ， 这 有 赖 于 管道 流量 的 增加 以 及 管道 温度 
的 升 高 ， 但 注意 到 管道 热 烟 损 在 数量 上 也 随 之 增 
加 ， 热 火 用 传递 效率 并 未 增加 。 此 外 ， 电 力 网 络 
有 籼 效 率 初 始 值 为 93.1 %， 下 降 稳定 至 92.6 %， 
减少 了 0.5 % 。 灯 合 元 件 烟 效 率 从 57.8 % 上 升 到 
60.5 %， 增 加 了 2.7% 。 


0.99 0,61 
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一 一 基态 网 络 : 
姓 0.97 一 电力 网 络 
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14 热力 / 电力 网 络 以 及 能 量 耦 合 元 件 中 的 籼 效 率 
Fig. 14 Exergy efficiency in thermal /electric power 


network and energy coupling elements 


分 析 可 知 ， 能 量 耦 合 元 件 烤 效率 不 到 61 
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%， 这 就 意味 着 流动 烟 价 至 少 比 电价 高 出 10.61 。 ”对 综合 能 源 系统 优化 运行 与 规划 以 及 综合 能 源 服 
倍 ， 当 能 量 看 合 元 件 效 率 提升 2.7 % 时 ， 热 烟 效 。” 务 等 具有 深远 意义 。 

率 并 没有 得 到 提升 ， 而 电能 的 品位 却 远 高 于 热能 。” ”参考 文献 
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